x = {'t-2';'t-1';'t';'t+1';'t+2';'≥ t+3'};
lwid = 2;

x2 = [-1:0.01:7];
x3 = [1;2;3;4;5;6];
y = zeros(length(x2),1);

C1 = [0.010632
0.013773
0.014574
0.024882
0.031225
0.021127];
 
tC1 = [0.87
1.24
1.25
1.97
2.42
1.82];

EC1 = 1.64*(C1./tC1);

figure(1)
hold on
fill([x3;flipud(x3)],[C1-EC1;flipud(C1+EC1)],[.8 .8 .8],'linestyle','none','FaceAlpha',0.3);
p1 = plot(x3,C1,'k','LineWidth',lwid);
plot(x2,y,'k')
set(gca,'FontSize',13)
ylabel('Coefficient','FontSize',19);
xlim([0.5 6.5])
ylim([-0.025 0.06])
xticks([1:1:6])
yticks([-0.02:0.01:0.06])
grid on
xticklabels(x);
hold off
print(gcf,'Figure5.tiff','-dtiff','-r300'); 